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Abstract 

The aim of this paper is to present an analytical calculation of the 
chemical potential of a Lennard Jones fluid. The integration range is 
divided into two regions. In the small distance region, which is r < a 
in the usual notation, the integration range had to be cut off in order 
to avoid the occurence of divergences. In the large distance region, the 
calculation is technically simpler. The calculation reported here will 
be useful in all kinds of studies concerning phase equilibrium in a 
LJ fluid. Interesting kinds of such systems are the giant planets and 
the icy satellites in various planetary systems, but also the (so far) 
hypothetical quark stars 
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1 Introduction 



The aim of this paper is to present a calculation of the chemical potential of 
a fluid consisting of neutral atoms or molecules. Interest in such systems has 
considerably increased towards the end of the last century, as a consequence 
of progress in planetary science. Until the end of October 2010, according 
to data at http://exoplane~teu| as much as 493 planets outside the Solar 
system have been detected. It has been shown that 423 of them have masses 
M < 5Mj where Mj is the mass of Jupiter. The major part (398) of the 
stars which have planets have masses equal to or smaller than the solar mass, 
and 151 planet has semi-major axis of the orbit between 1 and 3 astronomical 
units. 

Judging by experience from our Solar System, it is expectable that this 
interval of distances from a star corresponds to temperatures under which 
fluids consisting of neutral atoms and molecules can exist. It is known that 
giant planets have huge atmospheres and dense fluid interiors. Another class 
of planetologically interesting systems, in which fluids are important, are the 
icy satellites in our planetary system. For example, it has been concluded from 
data accumulated in the course of the Galileo mission, that jovian satellites 
Europa and possibly Callisto almost certainly have fluid oceans beneath their 
surfaces. 

Calculations to be discussed in this paper can find applications in the- 
oretical studies of quark stars. These are (so far hypothetical) phases of 
extremely dense matter, expected to occur in the interiors of neutron stars. 
In a recent study, aiming to constrain the parameters of solid quark matter 
by using data on the binary pulsar PSRJ1614 — 2230, the Lennard- Jones 
model was used to describe cold quark matter in quark stars [1] . It was 
shown there that if the number of quarks in one quark clusters is N q < 10 3 
there is enough parameter space for the existence of quark stars with masses 
higher than 2 solar masses. 

Modelling theoretically the internal structure of celestial objects ranging 
from the icy satellites and/or the giant planets to quark stars, demands the 
knowledge of the chemical potential of the fluid which they contain in their 
interiors. 

A necessary preparatory step in such a study must be the determination 
of the interparticle interaction potential. Obviouslyfor a fluid or any other 
kind of a system to be in equilibrium, the interparticle potential must be a 
combination of an attractive and a repulsive term. 
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It is known that between a pair of neutral atoms or molecules at a mutual 
distance larger than their respective dimensions there exists an attractive 
force - called the van der Waals (vdW) force (for example [2] or [3] ). The 
potential corresponding to the vdW force is proportional to r~ 6 ,where r is 
the interparticle distance. As shown by F.London, the physical origin of the 
vdW forces is the interaction of instantenous multipoles, while the repulsive 
contribution is of electrostatic origin. 

The vdW forces are anisotropic, which renders them additionally compli- 
cated [3] . However, their isotropic part is often approximated by the so called 
Lennard- Jones LJ potential. All the calculations in the following will deal 
with this particular model potential. The LJ model potential has the form 

(1) 

The symbol e denotes the depth of the potential well, while a is the diameter 
of the molecular "hard core". Obviously, lim r ^,ou{r) = oo. It can simply be 
shown that lim r ^ a u(r) = and that (du(r)/dr) = for r min = 2 1 / 6 <r. The 
depth of the potential well is u(r min ) = — e. 

An example of the LJ potential drawn for the particular case of CH^,with 
values of e and a from [5], is represented on figure 1. This particular molecule 
is interesting in two research fields: planetary science, because it is present 
in the atmospheres of the giant planets, but also in studies of the interstel- 
lar medium. On the figure,the distance is expressed in units of a and the 
potential divided by the Boltzmann constant ks is given in the units of K. 

2 The method of calculation of the chemical 
potential 

The chemical potential of a fluid (or any other system) can be calculated in 
two different ways: by using the general thermodynamical formalism,or by 
the general formalism of statistical mechanics. 

2.1 The thermodynamical formalism 

In this approach the calculation starts from the definition of the Gibbs po- 
tential: 

G = U - TS + PV (2) 
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Figure 1: The LJ potential for methane (CH4) 



where all the symbols have their standard meanings. Using the virial expan- 
sion, the pressure can be expressed as [4]: 

N 

P = P ID {\ + -B) (3) 

where N is the number of particles in the system and V the volume. The 
symbol B denotes the second virial coefficient, given by: 

1 r°° 

B = - (f - exp- uir)/T )dV (4) 
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and Pjd is the pressure of the ideal gas. The symbol u(r) denotes the inter- 
action potential. Inserting eq.(3) into eq.(2),it follows that 

G = U -TS + PV = G ID + NP ID B (5) 

The chemical potential is defined as \i = (dG/dN)p y T, which implies that 

^ = )p ' T + PidB + NB{? W )p > T (6) 

or 

dP 

fi = Vi D + Pi D B + NB(^) P , T (7) 



The equation of state of the ideal gas is PmV = NT which finally leads to 

fi = fi ID + 2P ID B (8) 

For the particular case of the LJ potential, it can be shown that the second 
virial coefficient is given by [5]: 

B(T) = -(b /2)fW^)(^ (9) 

where bo = 2na 3 /3 and T denotes the Gamma function. Inserting eq.(9) into 
eq(8) it follows that 
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h = hid- b oPlD E ^ r (^V i )(^)^ (io) 

n=0 U - 4 1 

which is the result for the chemical potential. Limitng the sum in this ex- 
pression to terms up to and including n = l,it follows that: 



A* - Hid + 2p ID b 



2.45083 - 1.8128(4) 1/2 



x(^) 1/4 (11) 



T 
T 

2.2 The formalism of statistical mechanics 

The chemical potential of a fluid is given by [6] : 

T ^- = Hp\') + -P- f 1 d 1 Tdr 
Kb-1 KbJ- Jo Jo 

xA7rr 2 u(r)g(r) (12) 

where 7 denotes the "coupling parameter" [6],p is the particle number density 
,u(r) the interaction potential and g(r) is the radial distribution function. 
The symbol h is the Planck constant divided by 27r,m is the particle mass 
and A is the thermal wavelength given by 

Expression (10) is valid under the condition pA 3 > 1, which leads to: 



The radial distribution function is a "bridge" relating macroscopic ther- 
modynamic properties with interparticle interactions in any kind of a sub- 
stance. In the theory of liquids, g(r) can be determined from first principles 
[6] just assuming a suitable form of the intermolecular potential [7]. In the 
following, the result for g(r) obtained in [7] will be used. Changing the 
variable from r to x = r/cr, and performing the integration over 7, it follows 
that ^ 

= ln( P A 3 ) + 47r<7 3 ^ jT dxx 2 u{x)g{x) (14) 

The domain of integration can be divided into two subdomains: x G [0, 1] 
and x G [1, 00], which means: 

poo 

I = a 3 dxx 2 u(x)g(x) = 
Jo 

pi roo 

a 3 [ / dxx 2 u(x)gi(x) + / dxx 2 u{x) x #2 (a?)] 

= a 3 x [h + h\ (15) 

This divergence of the LJ potential which occurs when x — > can be bypassed 
either by introducing a suitable change of the range of integration 16 [xq, 1] 
instead of x G [0, 1] with x 7^ 0,or by changing the form of the potential in 
the domain x G [0, 1]. For x G [0, 1] the function g(r) has the form 

gi(x) = sexp[— (mi + n) ] (16) 

and for x G [1, 00] the radial distribution function is 

g 2 (x) — 1 + \ exp[-(ax + b)] sin[(cx + d)] + 



x 2 



-^exp[-(0z + /i)]cos[(A;a; + Z)] (17) 



x 2 



where a,b,c,d,g,h,k,l,m,n and s are functions of pressure, temperature and 
density given in [7]. 

The appropriate boundary conditions, namely that the radial distribution 
function should tend to 1 in the limits of zero density and infinite distance, 
and the consequences of these conditions are also discussed there. As a 
consequence, the functions b, d, h and I are functions of density only, n is 
the function of temperature only and the other functions depend on the 
temperature and density [7]. 
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3 The calculation 



3.1 The case x G [0, 1] 

With the change of variables x = r/a,the LJ potential gets the form 

u(x) = 4e[x- 12 - x- 6 ] (18) 

Inserting Eqs (16) and (18) into the expression for Ji in eq.(15),it follows 
that 

00 ( — 1') /"111 



I x=xo 

\4J, 



(mx + npdx (19) 

Performing the integrations, after some algebra, it finally follows that 

8 27 
Ji -se x [1 + — m 4 - 

y 5 

27 o 9m 12 1 n 4 

— m 8 + .. + 5 + ..] 20 

10 54 2x1 2x o 

3.2 The case x G [l,oo] 

In this case, the calculation of the chemical potential is more straightforward. 
Inserting Eqs. (17) and (18) into the expression for I2 in eq.(15), and perform- 
ing the integration,gives the following approximate result for the integral I 2 : 

8e a 5 a 3 c 2 

J 2 = -— + 7recos[rf] cosh[6][— - — 

4 

+ 24 - ^cos[c]cos[rf] + ...] (21) 

3.3 The chemical potential 

According to Eqs. (14) and (15) the chemical potential is given by 

li = k B T ln(pA 3 ) + 4irpa 3 (h + h) (22) 

where the first terms of Ii and I 2 are given by Eqs. (20) and (21). 

Inserting Eqs. (20) and (21) into Eq.(22),one gets a simple analytical ap- 
proximation for the chemical potential of a LJ fluid. 
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4 Discussion and conclusions 



In this paper we have obtained an approximate analytical expression for 
the chemical potential of a Lennard Jones fluid. Two ways in which such an 
expression can be obtained have been presented, and both of these approaches 
has been applied. 

The approach based on the general thermodynamic formalism gives a re- 
sult, expressed as eq.(10), which is both mathematically and physically sim- 
pler. It contains just two variables which characterize the material under 
consideration - these are a - the diameter of the molecular "hard core", and 
e - the depth of the potential well. Note that the chemical potential obtained 
in this way for a certain value of the ratio e/T reduces to the value liid- 

The formalism of statistical mechanics is both mathematically and phys- 
ically more complex. The general conclusion is that the chemical potential 
depends on the thermodynamic parameters of the fluid through the func- 
tions a-s, which are in turn functions of the pressure and/or density and/or 
temperature [7] ,but also on the interaction parameters. The approximate 
expression for the chemical potential of a LJ fluid is: 

39 97rr? 4 
Li = k B Tln(p\ 3 ) + ^npeah[l + 

9 5 

+.. + + .. + — cosfdl cosh[&l( — — ) 

2x 9 s L JV 320 32 ; 

+..] (23) 

All the symbols in this expression have their standard meanings, or were in- 
troduced in (Morsali et al.,2005). Mathematically, the symbol xq denotes 
the cut off radius of the LJ potential introduced in the calculations in order 
to avoid the occurence of divergences. Physicallythis quantity represents 
the interparticle distance at which pressure ionisation occurs. Qualitatively 
speaking, pressure excitation and/or ionisation occur because electronic en- 
ergies change under the influence of the external pressure field. For details 
about this process see,for example, [8]. 

The calculation presented in this paper was motivated by recent advances 
in planetary science. As a consequence of numerous discoveries of giant 
exoplanets, modellisation of their internal structure has regained importance. 
These planets consist mostly of fluids, and accordingly an obvious need for 
theoretically "preparing the ground" for the modellisation of their interiors 
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has occured. Studies of phase equilibrium and phase transitions demand an 
explicit knowledge of the chemical potential. Some preliminary results in 
that direction have recently been obtained [9] in the limit of small density 
and without taking into account the chemical potential. Another interesting 
problem,which becomes accessible for study with the results obtained in this 
paper is the behaviour of the chemical potential of a LJ fluid with changes 
of its thermodynamical parameters. Some aspects of both of these problems 
will be discussed in future work. 
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